Quantum Hall Effect and Semimetallic Behavior of Dual-Gated ABA-Stacked Trilayer Graphene 
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The electronic structure of multilayer graphenes depends strongly on the number of layers as well as the 
stacking order. Here we explore the electronic transport of purely ABA-stacked trilayer graphenes in a dual- 
gated field-effect device configuration. We find that both the zero-magnetic-field transport and the quantum Hall 
effect at high magnetic fields are distinctly different from the monolayer and bilayer graphenes, and that they 
show electron-hole asymmetries that are strongly suggestive of a semimetallic band overlap. When the ABA 
trilayers are subjected to an electric field perpendicular to the sheet, Landau level splittings due to a lifting of 
the valley degeneracy are clearly observed. 

PACS numbers: 73.22.Pr,73.20.At,73.43.-f,73.50.Dn 



Single and bilayer graphenes are well known to have low- 
energy electronic dispersions that differ dramatically from the 
parent compound, graphite fl]-[3]. With just one more layer, 
two stable forms of trilayer graphene are known to exist hav- 
ing either ABA (Bernal)- or ABC (rhombohedral)-stacking, 
which exhibit differing electronic structures. In particular, 
ABA-stacked trilayer graphene is predicted to be a semimetal 
and is thus the first- and thinnest- multilayer graphene to have 
a band structure resembling that of bulk graphite |4}j9]]. Not 
only do trilayer graphenes comprise a new class of materials 
intermediate between graphite and its fundamental building 
block of monolayer graphene, but studies of the electronic 
transport of ABA graphenes may also shed light on trans- 
port anomalies observed in bulk graphite and other semimet- 
als whose properties cannot be as easily varied by externally 
applied gate voltages ifTOlfTTl . 

Trilayer graphene has received increasing attention of late, 
due to interest in the unusual electronic structures of ABA- 
and ABC-stacked graphene, and also to advances in identi- 
fying the nature of the stacking order of trilayer flakes. The 
two crystal structures differ in that for ABA graphene the top 
and bottom sheets are aligned while the middle sheet is shifted 
by one bond length along the C-C bond direction (see Fig. 1 
(a) ); while for ABC graphene the top layer is also shifted, 
by two bond lengths relative to the bottom layer. Measure- 
ments of the electronic transport in trilayers have been previ- 
ously reported |[T2l . and the nature of the stacking has been 
inferred by comparing the observed signatures of the quan- 
tum Hall effect (QHE) to those expected for ABA |fl3] and 
ABC 1 14] graphenes. However, the recent development of 
optical methods to clearly distinguish between regions of dif- 
fering stacking order has lead to investigations of the elec- 
tronic properties of trilayers known to be chiefly of one type 
|[T5HT9l . In this work, we study large area dual-gated trilayer 
graphene samples known to be of ABA stacking and present a 
unified view of the electronic transport of this system that has 
not been available to date. We find that the zero-, low-, and 
high-magnetic-field quantum Hall transport of ABA-stacked 
trilayer graphene samples can be understood in the context 
of the underlying semimetallic band structure that arises from 
the mirror-symmetric ABA crystal lattice. 



On average, trilayer graphene contains regions of ABA- and 
ABC-stacking in an 85:15 ratio. The two forms pose a prob- 
lem for measurements of the QHE where the nature of the 
edge states may change, depending on the stacking type, to 
such an extent that even a small region of differing stacking 
order can affect the interpretation of transport data. Therefore 
in this work only those graphene flakes identified as being pre- 
dominantly composed of ABA-stacked trilayers by Fourier- 
transform infrared (IR) spectroscopy are utilized. Moreover as 
for bilayer graphene, in multilayer graphenes the full range of 
electronic transport cannot be accessed by a single-side gated 
device, since varying the voltage on just one gate simultane- 
ously changes both the Fermi level (or carrier density n) and 
the drop in electric potential (or electric displacement field, 
D) between the layers [20-22]. This problem is remedied by 
the use of a second (top) gate which allows for independent 
control over the values of n and D. In ABA-stacked trilayers 
we will see it is crucial to distinguish between effects that de- 
pend on only one or the other parameter in order to understand 
the observed electronic transport. 

In Fig. 1 (b) an optical microscope image of a typical device 
with an area of 3000 /im 2 is shown. To determine the layer 
stacking, each flake is studied via Fourier-transform IR spec- 
troscopy prior to any device fabrication. The IR reflectance 
at room temperature, i? t n, averaged over the entirety of the 
graphene sheet embedded in this device and normalized by 
the reflectance of the nearby substrate, i? su b, is shown in Fig. 
1 (c). Two features stand out: A broad minimum centered 
at 6000 cm -1 and a smaller dip near 4000 cm -1 . While the 
former is due to interference effects in the 285 nm-thick SiC<2 
layer 11231 . the latter feature arises from transitions to (from) 
the split-off conduction (valence) band. The relevant transi- 
tions when the system is lightly doped are shown schemati- 
cally in the inset to Fig. 1 (c). These are similar to the IR 
absorption in bilayer graphene although higher in energy by a 
factor of v2 124W261 . Along with the absence of any IR sig- 
nal due to ABC graphene which would appear at 2500-3000 
cm -1 , and with independent confirmation via Raman spec- 
troscopy measurements, the 4000 cm -1 dip clearly identifies 
the flake as ABA-stacked trilayer graphene lfT5l[T6l . In partic- 
ular, analysis of the reflectance data allows an estimate for the 
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FIG. 1. (Color) (a) Crystal structure of ABA (Bernal)- stacked trilayer graphene, with arrows indicating the tight-binding hopping parameters 
of the Slonczewski-Weiss-McClure model |8 9|. (b) Optical image of a 3000 /im 2 dual-gated ABA trilayer graphene device, (c) Infrared 
reflectance, R tI [, of the trilayer flake shown in (b) normalized to the substrate reflectance, -R sll b- The dip near 4000 cm -1 is due to transitions 
to or from the split off bands, as shown schematically in the lower left inset for the case of light electron doping. The larger dip at 6000 cm -1 
is due to interference in the thin Si02 substrate 1231 . 



content of the ABC-stacked trilayer graphene in this sample 
at less than 3%. 

Graphene-on-SiC»2 samples found to have ABA stacking 
are fabricated into dual-gated field-effect devices via stan- 
dard electron beam lithography and thin film deposition tech- 
niques. The top gate dielectric is composed of polymethyl 
methacrylate (PMMA) 950 K A4 diluted with methoxyben- 
zene (anisole), and exposed to a 25 kV electron beam at 
a dosage of 24, 000 /xC/cm 2 to create a hardened dielec- 
tric layer. Applying voltages Vb and V t to the back and 
top gates leads to changes in both n and D, as given by 
n = a (AV h + pAV t ) and I? = (D b -D t )/2 = (e h AV h /d h - 
e t AV t /d t )/2, with AV- = (Vi - V i0 ). Here the d[ and e ; are 
the thicknesses and dielectric constants of the insulating lay- 
ers, respectively, and the Vjo are the offset voltages required to 
reach n = D = due to extrinsic doping. Both the density 
and the gating efficiency, a = 7.5 x 10 10 cm _2 V _1 , are cal- 
ibrated by oscillations in the magnetoresistance, R xx , at high 
magnetic fields. The ratio of the two gate capacitances, j3, 
is determined from the slope, to, of the high-resistance ridge 
in the zero-field resistivity measurements in Fig. 2 (a) by 
—to = 1//3 = d t e\, / (d^Ct) . The precise location of D = 
is determined via features in the QHE, as described below. 
All electronic transport measurements are made at T = 0.3 K 
using standard low-frequency lock-in techniques. 

In Fig. 2 we show the sheet resistivity, p, at zero magnetic 
field, B, as a function of the gate voltages V t and Vj,. Two 
arrows superimposed on the data define the axes of increasing 
n and D induced by the two gates; they cross at a saddle point 
in the resistivity which is independently identified as where 
n = D = via measurements at high magnetic field. This 
saddle point bears a superficial resemblance to the gate volt- 



age dependence of the resistivity of bilayer graphene l27ll28l . 
Overall, the transport here is outlined via profile cuts along 
the D and n axes, plotted in Fig. 2 (b) as the conductivity, 
<7, against n (top axis) or D (bottom axis) respectively. Data 
from two samples are shown, with the solid lines taken from 
Fig. 2 (a), and the dotted lines from a second sample in which 
the saddle point was further displaced from Vb = V t = 
due to extrinsic doping, thereby allowing a greater range of 
D to be accessed. As for monolayer and bilayer graphenes 
on SiC<2, a oc n away from charge neutrality. At the saddle 



point, <7 n 



= 7- 



/h and is 2 to 3 times greater than for 



typical bilayer graphene on SiC»2, perhaps because the higher 
density of states in trilayer graphene leads to more effective 
screening of scattering sources [29-31 1. The change in a with 
increasing D, however, is the most unique feature. Initially, a 
decreases by a factor of 2 or 3 until D/e a ss 0.4 V/nm, beyond 
which it is flat until, for D/eo > 1.0 V/nm, it begins to slowly 
rise again. This new and unusual behavior- observed in mul- 
tiple samples- strongly suggests that, unlike bilayers, no gap 
opens in ABA trilayers at high D. This behavior is similar 
to predictions for the conductivity of ABA trilayers subject to 
an electric potential imbalance between the top and bottom 
graphene sheets, in which a(n = 0) is generically expected 
to decrease sharply and then rise slowly as D increases from 
zero ID. 

We note that earlier work on dual-gated trilayer devices 
found entirely different behavior for the resistivity, show- 
ing an apparent global maximum at the point identified as 
n = D = H2M19). However, the location of the D = 
point was not independently determined in those studies. 

At high magnetic fields, the QHE for ABA-stacked tri- 
layer graphene exhibits a number of properties that distin- 
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FIG. 2. (Color) (a) ABA trilayer sheet resistivity, p, in units of kf2 per square, vs. the top and back gate voltages, V t and Vb. The arrows at the 
saddle point indicate the axes of increasing carrier density, n, and electric displacement field, D. (b) Profile cuts along the n and D axes are 
plotted as the conductivity, a, vs. n or D. The solid lines are from Fig. 2 (a), while the dashed lines are from a second sample in which the 
saddle point is offset from Vb = Vt = due to extrinsic doping, so that a larger range of D is accessible. 



guish it from other graphenes. In Fig. 3 (a), the inverse of 
the measured Hall resistance, R~y, is plotted for B = 14 T. 
Several QH plateaus are clearly visible as wide stripes run- 
ning from the upper left to lower right. Contour lines drawn 
at half-integer values of the filling factor, v = nh/eB = 
... — 1.5, —0.5, +0.5, +1.5..., serve to emphasize the boundary 
of each plateau so that, for instance, the —6 plateau is bounded 
by contours at v = —6.5 and —5.5. The interesting features 
of this data fall into three groups: First, the most obvious 
plateaus, running diagonally in parallel to the D axis, com- 
prise the sequence of quantum numbers —14, —10, —6, +6, 
and +10 that are labeled in the figure. Second and most un- 
usual, several distinct plateaus with quantum numbers ±2 and 
±4 develop near the upper left and lower right corners of Fig. 
3 (a), at low n and high values of D. Finally, a small region 
at the center of the plot is found to be well quantized with a 
conductance of —2e 2 /h. For all of these, the precise quanti- 
zation of the plateaus is highlighted by cuts made along the 
n (D — 0) axis, as well as at the boundaries of Fig. 3 (a) at 
constant V t for varying Vb, and vice versa, all shown in Figs. 
3 (b)-(d), respectively. 

Several of the QH features can be explained with recourse 
to the lowest-order model of the ABA trilayer band structure. 
In a tight-binding calculation when only terms describing the 
in-plane nearest-neighbor hopping and the largest interlayer 
hopping- 70 and ji, respectively, in Fig. 1 (a)- are included, 
the resulting band structure arising at the K and K' points 
in momentum space is depicted in Fig. 4 (a). It resembles 
a "1+2" superposition of linear and hyperbolic bands similar 
to those of monolayer and bilayer graphene |5j [32) . On the 



right of Fig. 4 (a) the corresponding zero-energy Landau lev- 
els (LLs) are schematically plotted as a function of position, 
so the electron or hole levels diverge upward or downward to- 
ward the sample edge. In monolayer (bilayer) graphene, one 
(two) fourfold degenerate LL(s), reflecting the spin and valley 
degrees of freedom, are found at zero energy. Thus, in this 
1+2 model, three such zero-energy levels appear and give rise 
to an overall 12-fold degeneracy, yielding quantum numbers 
+6 or —6 when the Fermi level lies where indicated by the 
dashed lines. Consequently, accounting for additional higher- 
energy LLs will, barring accidental degeneracies, lead to the 
sequence ... — 14, —10, —6, +6, +10, +14... which is identi- 
cal to that observed along the line where D = except for the 
presence of the —2 plateau. 

The appearance of the ±2 and ±4 quantized plateaus in 
the upper left and lower right corners of Fig. 3 (a), high- 
lighted in Fig. 3 (c) and (d), can also be understood in a low- 
est order model of ABA trilayers. In the theoretical picture 
of Ref. [32 1, the ABA trilayer Hamiltonian may be decom- 
posed into one monolayerlike and one bilayerlike system, and 
since each of these alone has inversion symmetry, an effective 
inversion symmetry is recovered for the ABA trilayer. This 
guarantees the degeneracy of states at the K and K' valleys, 
despite the lack of a true lattice-inversion symmetry. How- 
ever, this model does not hold in the presence of a potential 
drop between the graphene layers due, e.g., to the external D 
field, and the valley degeneracy is thus lifted. In particular, 
when D / 0, the three fourfold degenerate LLs at zero en- 
ergy split, so that one twofold degenerate level (for electron 
spin) from each valley remains at zero energy, while the re- 
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FIG. 3. (Color) (a) The ABA trilayer quantum Hall effect, plotted as the inverse Hall resistance, R~} , vs. VE, and V t at B — 14 T. Contour 
lines are drawn at filling factors v = ... — 1.5, —0.5, +0.5, +1.5.... The quantum numbers of several plateaus are indicated in white. Several 
profile cuts through the data to highlight the various quantized plateaus are shown in (b), (c), and (d). The cut in (b) is taken along D = 
following the dashed line in (a). The cuts for (c), for constant top gate voltages, are taken along the top and bottom edges of (a). The cut in (d) 
is taken along the right edge of (a). 



maining levels disperse with a dependence on D that differs 
in each valley. Therefore, as D increases from zero, these di- 
verging LLs will initially give rise to new quantized plateaus 
in a 4,2,-2,-4 sequence ll32l[33l that is strikingly similar to the 
pattern of the QH plateaus observed in Fig. 3. With further 
increase of D, additional plateaus are expected at even val- 
ues of the filling factor, and indeed in a second sample the 
v = +8 and +12 plateaus have also been seen. Since only 
the strength, and not the sign, of the applied field is important, 
the 4,2,-2,-4 structure appears symmetrically for ±D. Thus, 
this feature can be used as an independent method to locate 
the line where D = 0, which is found to pass through the sad- 
dle point of Fig. 2 (a), justifying the previous identification of 



the saddle point as where n = D = 0. Finally, we note that 
the QHE in ABC trilayers for devices of similar quality is ex- 
pected to be distinctly different, with a filling factor sequence 
v= ... - 6. 0. +6... for any value of D l34t 

While the lowest-order model of ABA trilayers successfully 
explains several of the QH plateaus in Fig. 3 (a), it offers no 
insight into the presence of the robust —2 plateau near the cen- 
ter of the plot, nor to the corresponding asymmetry reflected 
in the lack of a plateau at v = +2. Clearly, the v = — 2 
plateau must arise from a splitting of the three zero-energy 
LLs. While in principle many-body interactions can be re- 
sponsible, in these low-mobility samples (/i 4000 cm 2 /Vs) 
disorder is expected to overwhelm such effects. Therefore, 
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FIG. 4. (Color) (a) The band structure of the "1+2" model of ABA 
trilayer graphene (left) and the associated lowest LLs (right). QH 
plateaus will occur when the Fermi level lies at the dashed lines, 
with quantum number +6 (—6) when all levels are full (empty), (b) 
Band structure calculated to higher order in the Slonczewski-Weiss- 
McClure model (see text). Two bands overlap at zero energy, creat- 
ing a semimetal. To the right, the corresponding lowest LLs allow 
for additional QH plateaus at —2, +2, and +4. Level energies are 
plotted against position, diverging upward (downward) for electron 
(hole) states near the sample boundary. All LLs are twofold degener- 
ate for the electron spin. Closely spaced pairs only mark the number 
of spin-degenerate levels, not additional splittings. 



we seek an explanation for the —2 plateau in terms of a more 
accurate model of the trilayer band structure ll3~5ll . 

In Fig. 4 (b), the ABA band structure is shown, calculated 
to higher order by utilizing the full set of tight-binding param- 
eters in the Slonczewski-Weiss-McClure model for graphite 
1 1 1 with values identical to those used in Ref. JS): vq(<x 70) = 
1.0 x 10 6 m/s, 71 = 0.4 eV, 72 = -O.O571, v 3 (<x 73) = 
O.lurj, w 4 (oc 74) = 0.014v , 75 = O.I71, and S = O.I2571, as 
well as D = 0. In this picture, the monolayerlike and bilayer- 
like bands become gapped as well as offset from zero energy, 
leading to the overlap of one approximately linear hole band 
with a nearly parabolic electron band. The region of overlap 



delimits the semimetal: At charge neutrality, both bands are 
partially occupied and contribute to transport. 

The lowest energy LLs for the higher order calculation are 
plotted on the right of Fig. 4 (b) as twofold degenerate levels 
(for electron spin), with one (two) arising at each extrema of 
the monolayerlike (bilayerlike) bands. Some surprising fea- 
tures appear: New gaps develop at v = —2, +2, and +4, and 
a LL crossing naturally arises as the levels diverge at the sam- 
ple edge. Thus, the higher order model can account for the ob- 
servation of a plateau at v = —2. But why then are plateaus 
absent at v = +2 and +4? In fact, the relative sizes of the 
gaps shown at v = —2, +2, and +4 in Fig. 4 (b) may not be 
correct: The gaps shown here depend sensitively on the val- 
ues of the tight-binding parameters used in the higher-order 
calculation, and values of these are not well known for tri- 
layer graphene. The data in Fig. 3 would suggest that the 
v = — 2 gap is larger than those at v = +2 or +4, which 
may be washed out by disorder. However we note that a dip 
appears in f?~* near u=lx 10 12 cm~ 2 in Fig. 3 (b), pos- 
sibly signaling an incipient +2 or +4 plateau. Figure 4 (b) 
further illustrates an intriguing feature of the v = +2 state, 
due to the LL crossing of electron and hole bands which carry 
counter-propagating edge modes. While it is unclear whether 
this should affect the quantization, this feature may lead to 
new physics in future high mobility devices I36ll37l . 

As the LL sequence and electron-hole asymmetry in the 
center of Fig. 3 (a) are both consistent with a semimetallic 
band overlap in ABA trilayers, it is natural to assume that 
low-field measurements will also show semimetallic trans- 
port. To address this assumption, the low-field Hall coeffi- 
cient, i?H = dR X y/dB, and sheet resistivity, p, are compared 
to a control sample of monolayer graphene. The results are 
shown in Fig. 5 (a) and (b) for trilayer and monolayer sam- 
ples, respectively, where the measured R H and p are plotted 
against carrier density. The trilayer data are acquired along 
D = 0, which is possible only in dual-gated samples. The 
evolution of J?h is quite different for the two systems. The 
monolayer data are perfectly symmetric about n = 0; due to 
disorder, they do not diverge but rather smoothly change sign 
over a narrow region at low densities. In contrast, the trilayer 
data are strongly asymmetric, having a large positive peak at 
n = and a zero crossing at slightly more than n = 1 x 10 12 
cm~ 2 , followed by a weak minimum and slow decrease back 
toward zero. Meanwhile, in ABA trilayers, a p xx ~ B 2 
behavior is observed, characteristic of two-band conduction. 
This observation, along with analysis of the clear two-band 
signatures in the Shubnikov-de Haas oscillations, will be dis- 
cussed in future work. 

The different behavior of R H can be qualitatively under- 
stood in terms of semimetallic transport. As expected for con- 
duction arising from a single band, the monolayer data are 
well described by Ru = —1/ne except at the lowest densi- 
ties, which are modeled by convolving Ru with a Gaussian 
whose variance, 5n, represents a spread of densities due, e.g., 
to electron-hole puddles IF2T1I381 . While this model certainly 
oversimplifies the role of disorder, a reasonable fit to the data 
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FIG. 5. (Color) (a) The low field Hall coefficient, Rh, and sheet resistivity, p, in kQ per square, measured in ABA trilayer graphene as a 
function of n, for D — 0. Uncertainty in the measurements is smaller than the symbol size, (b) Rh and p measured in monolayer graphene. 
The dotted line is Rh = —1/ne, averaged over disorder (see text), (c) Simplified ABA band structure, (d) Rh vs. n calculated for the band 
structure in (c), for mobility ratios pi/p p = 3.3 (solid blue), 1 (dashed red) and 0.3 (dotted green) (see text). 



can be found for 5n = 0.8 x 10 11 cm~ 2 (dashed line in Fig. 5 
(b) ). In contrast, the transport properties are qualitatively dif- 
ferent when two bands take part in the conduction, as occurs 
in ABA trilayers. In particular, when electron and hole bands 
overlap, the low-field Hall coefficient becomes B9l 

1 p-n (^e/Mh) 2 m 

-n-H = 2 ' ' ) 

e (p + rifle/ Hh) 

where p and n are the hole and electron densities and fte/fih 
is the ratio of their mobilities. Sharply differing from the 
behavior for a single band, this two-band form of i?n will 
in general cross zero away from charge neutrality (p = n); 
and at the zero-crossing, the ratio of the densities is given by 
p/n = (fijfi h ) 2 . 

In Fig. 4 (b), the number of conducting bands- and as- 
sociated mobilities- changes several times as the Fermi level 



is varied. To model R^, the simplified ABA band structure 
shown in Fig. 5 (c) is adopted. Depending on the Fermi 
level position, the relevant one- or two-band expression is 
employed, and the result is smoothed with a Gaussian using 
6n = 3 x 10 11 cm" 2 . Obviously, these results will vary with 
the choice of mobilities and band offsets, but overall the be- 
havior of this model can be summarized by the curves in Fig. 
5 (d). There, Ru vs. n is calculated for three ratios of the band 
mobilities, fi\/ fi p — 3.3, 1, and 0.3, where fi\ and fi p are mo- 
bilities of the linear and parabolic bands, irrespective of the 
sign of the carriers. Interestingly, the data are best approxi- 
mated when fi\ > /i p , which gives the only curve (blue line 
in Fig. 5 (d) ) that captures the pronounced asymmetry of the 
data and crosses zero at a positive carrier density. Of course, 
the many choices of band-edge offsets and curvatures, com- 
bined with the likely density-dependent mobilities, all prohibit 



7 



using this model to make accurate fits to the data. In particu- 
lar, substrate-induced potential fluctuations are ubiquitous in 
SiC>2 -based devices and likely to contribute to a smearing of 
transport effects over a range of densities |3]|40]. Nonethe- 
less, the simple semimetallic model of Fig. 5 (c) does capture 
the essence of the data. These results, in concert with the 
clear asymmetry of the quantum Hall effect, strongly suggest 
that ABA trilayer is indeed a semimetal with an unusual LL 
structure at low energies. 
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